T-duality in supersymmetric theory of disordered quantum systems 
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A new super-symmetric representation for quantum disordered systems is derived. This represen- 
tation is exact and is dual to that of the nonlinear sigma-model. The new formalism is tested by 
calculating the distribution of wave function amplitudes in the Id Anderson model. The deviation 
from the distribution found for a thick wire is detected near the band center E = 0. 
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Introduction. — The rigorous theory of quantum dis- 
ordered systems is based almost entirely on the field- 
theoretical formulation of the nonlinear sigma-model (NL 
<tM) 0,0. The most successful variant of this theory, 
the so called supersymmetric (SUSY) sigma-model 
allows not only for a perturbative treatment of weak lo- 
calization and mesoscopic phenomena in the delocalized 
(metallic) phase which was previously done by the impu- 
rity diagrammatic technique |^| but also applies to sys- 
tems in the localized (insulator) phase provided that the 
localization radius is sufficiently large. A celebrated ex- 
ample is the problem of localization in a thick disordered 
wire 0, HI which is equivalent to the Gaussian ensem- 
ble of random banded matrices with a large bandwidth 
&> 1. 

However, despite of the great success of the nonlinear 
SUSY sigma-model over the last two decades Q this ap- 
proach is limited. The limitation comes from the saddle- 
point approximation that locks the norm of the super- 
matrix field Q and thus leads to a geometric constraint 
Q 2 = 1. It is precisely the validity of this saddle-point 
approximation that requires a quantum system to be in 
the metallic phase (at least for a sufficiently small length 
scale), or the banded random matrices to have a large 
bandwidth. 

The best known example for which the NLcrM is in- 
sufficient is the strictly one-dimensional disordered chain 
described (in the position representation) by the three- 
diagonal matrix Hamiltonian fljj = Ji 5j^±± + Ei §ij. The 
case Ji — 1, and £j being a Gaussian random variable 
with zero mean and the variance w = (ef) is known as 
the Id Anderson model. Its continuous counterpart for a 
weak disorder w <C 1 was studied in detail in the 70-th 
0, [M with the classic results on the frequency depen- 
dent conductivity [tI|10|| and dc conductance distribution 
H . The special case where = and Ji is a Gaussian 
random variable has been studied in the pioneer work by 
Dyson [ijj and attracted a lot of interest in recent years 
in connection with the new symmetry classes fl2|. Yet 
another broad class of systems which cannot be studied 
by the NLcrM are systems described by the almost diag- 
onal random matrix Hamiltonian |l.'{l| with the variance 
of the diagonal entries (Hfj) much larger than that for 



the off-diagonal ones (\Hi^j\ 2 ). Some of such systems 
with (\Hi^j\ 2 ) ~ \i — j\~ 2 attracted recently much atten- 
tion in connection with the multifractality of the critical 
eigenstates [lHHHlH. 

The main goal of this work is to fill in the gap in our 
technical ability to apply the SUSY method to quantum 
disordered systems. Below we present a field-theoretical 
description that is free of of the saddle-point approxima- 
tion and thus applies to all the above cases which are 
not amenable to the conventional NLcrM. This descrip- 
tion keeps all the Grassmann variables and thus is dif- 
ferent from another suggestion, the su per - opsonization, 
which was recently proposed in Ref. |l6| as an exact 
method instead of the NLcrM but is suffering from un- 
certainties and incomplete definition of the target space. 
The method we are proposing is similar in spirit to that of 
Ref.0 but has an advantage to be as close as possible to 
the standard SUSY NLcrM. In particular, the new super- 
matrix field Q can be diagonalized by a pseudo-unitary 
transformation with separated Grassmann and commut- 
ing variables Q similar to the one used in the NLcrM. 
This allows to drastically simplify calculations with the 
final result depending on few commuting variables only, 
as for the conventional NLerM. 

The main property of the new description is that it 
is dual to the old one: the coupling matrix of the new 
theory is just the inverse of the coupling matrix of the 
NLcrM. However, the space of the matrices Q of the new 
theory (the Q target space) is different from the Q tar- 
get space of the NLcrM. So, we are dealing not with the 
conventional duality but with a kind of T-duality [l7^ 
remotely resembling the situation in the string theory. 

As an illustration of the new formalism at work we 
consider the problem of the eigenfunction statistics for a 
one-dimensional Anderson model. We derive the trans- 
fer matrix integral recursive equation following the same 
rout as the one used in the framework of the SUSY non- 
linear sigma-model for quasi-ld disordered wires 0] or 
banded random matrices 6]. For weak disorder this in- 
tegral equation reduces to the the differential equation of 
exactly the same form as for a thick disordered wire. This 
fact (noticed earlier in ) implies that despite the ab- 
sence of the diffusive regime in a strictly Id chain, the 
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smooth envelope of a wavefunction at scales much larger 
than the wavelength and the wire cross-section has ex- 
actly the same statistics in a thick disordered wire and 
in a strictly Id disordered chain. 

The dual representation. — For concreteness we con- 
sider the matrix Hamiltonian of the unitary class H 
which entries Hij are random independent Gaussian vari- 
ables characterizing by the variance matrix gy = (|-ffy | 2 ) 
and the mean values (Ha) = . The standard 
SUSY field-theoretical description 3] aimed to com- 
puting the retarded/advanced (R/A) Green's functions 
[E± — H] _1 begins with introducing the W functional 
S = Eij *i (E± Sij - Hij) where E± = E ± Cj/2 
and Q = lu + 2ir\ with r\ — > +0. The ket -supervector 

= {Sr, Sa, M-Ri I jL a) T contains both commuting com- 
plex variables Sri a and Grassmann anti-commuting vari- 
ables Hr/a- For convergence of the Gaussian functional 
integrals containing e lS ^ the bra -supervector ^/ is given 
by * = *tA where A = diag{l, -1, 1, -1}. The SUSY 
trick allows to do averaging over disorder prior to do- 
ing the functional integral over 'J'. Thus we obtain 
(e zS [*]) = e~ F , where 

F = \ E &i ^AQiQA - y £ Str t A ^ C 1 ) 

- iE£Str[&] + F Wn . 

Here F fcm = i £y *< Hf ^ and 

Q = $® ^ = 4r® *t a. (2) 

So far we did not deviate from the standard rout ex- 
cept for introducing the new notations Eq.©. How- 
ever, the integration measure is still unchanged (i/x^] = 
dS roIS AdS' R dS'\d^RdpLAd^ R d^jL A . 

Now we do the crucial step. Instead of performing 
the Hubbard-Stratonovich transformation which decou- 
ples the quartic in ^ term in Eq. © by the supermatrix Q 
and applying the subsequent saddle point approximation 
which leads to the sigma-model constraint Q 2 = 1 we 
change the variables so that in the new variables the su- 
permatrix Q is parameterized similar to the supermatrix 
Q of the sigma-model: 



~>R/A 



±iJ^~2e ±i */ 2+in (l-lxR/AXn/A) (3) 



Vr/a = ±wAi/ 2 e 



XR/A, 



where A1/2 > 0, < ip < 2tt, < f2 < n and Xr/a, X*r/a 
are the new Grassmann variables. 

One can check that in new variables the supermatrix 
Q defined in Eq. (JSj can be represented in the form where 
Grassmann and commuting variables are separated by 
factorization: 
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(4) 



The factorization Eq.0} is also one of the basic proper- 
ties of the Efetov's parameterization in the sigma-model 
0. Moreover, the corresponding matrices Ur/a are es- 
sentially the same: 



ur/a 



jXr/aXr/A 
Xr/a 1 



~X*R/A \ 



(5) 



2Xr/aXR/Aj 



BF 



The difference is in the form of the matrix D. Like in 
the case of Efetov's parameterization it is diagonal in the 
boson-fermion (BF) space. However, in contrast to the 
sigma-model only the bosonic sector is nonzero: 



D BB = 



Ai 



-yfXTMe-^ -A 2 



D FF = 0. 



RA 



(6) 

One can see that the matrix Ebb can be diagonal- 
ized by the same pseudo-unitary transformation as the 
matrix Ebb of the sigma-model: Ebb = T BqT^ 1 , 
where T G 17(1,1). However, the diagonal matrix 
B = |Ai - A 2 |diag(0(Ai - A 2 ),0(A 2 - Ai))** is differ- 
ent from the corresponding diagonal matrix Abb of the 
sigma-model and contains a new degree of freedom, the 
difference Ai — A 2 € R. This means that the full sym- 



metry of the target space is jjtj 



u(i,i) 



U(l,l) 



U(2) 



(l)xU(l) 



rather than 



as for the sigma-model. 



U(l)xU(l) ~ U(l)xU(l) 

To complete the procedure of changing variables from 
f to Q we have to compute the Jacobian of the transfor- 
mation Eq.© and rewrite Fki n in Eq.© in terms of Q. 
The first task is straightforward and the result is: 



d\i = dQdfl = dxRdxRdxAdx*Adipd\id\2 



dn 

AiA 2 



(7) 



To accomplish the second one we note that = <pi e 



and 



where tpt = ^ 



i\n i= 0- 



Then each 



link {ij} enters in the last term of Eq.© as iVij = 
i(Pi 



-iQij 



Pji e iS! «), where p i3 = (p l (p j H^°\ and fl 



Qi — Qj . Then assuming that the observable of interest 
is gauge-invariant and thus independent of fl and inte- 
grating e~ l ^i3 Vi i over all fly £ [— ir,ir] one obtains for 
the last term in Eq.©: 



j>i 



In i 2tt Jo 2 



H (0) | 2 



Str[QiQj 



(8) 



Note, however, that in order to arrive at Eq.© inte- 
gration over all VLij should be independent. This is the 

case when the links provided by form a lattice with 
no loops. The simplest example considered below is a 
one-dimensional lattice with only nearest neighbor con- 



nections H, 



(0) 



i,i±l- 



The new supermatrix field theory given by Eqs.©,© 
together with the parameterization of the supermatrix 
Q, Eqs.QJ-© is the main result of the paper. It has to 
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be compared with the standard supermatrix sigma-model 
In the case (fly) = and gy = g(| i — j |) the action 
reads lol: 



(iri> A ) 



i 

(9) 



where v{E) = j^-t/SAq - E 2 is the mean density of 
states and A = £\ gy . Comparing Eq. @ with Eq. 
one notices that the coupling matrix g in Eq.Q is re- 
placed in Eq.© by the inverse coupling matrix g -1 . 
That is why the new representation is dual to the stan- 
dard sigma-model representation. However the target 
spaces of the supermatrices Q and Q in Eq. and Eq. © 
are different. The Q-field is constrained: Q 2 — 1 and 
Str Q = 0. That is why the term E Str Q cannot arise and 
the energy E enters as an inessential parameter. There 
is no such constrains for the supermatrix Q in Eq.l|T|) 
and the terms EStrQ and giiStrQ 2 are important. In 
particular this leads to peculiar prop erties of the Id disor- 
dered chain at the point E = [ljj, |2J| , especially in the 
chiral case where diagonal matrix elements (on-site 
energies) Ha do not fluctuate. 

Id Anderson model. — We illustrate the new method 
by considering the statistics of an eigenfunction ip n (i) 
corresponding to the eigenvalue E n for the Id Anderson 
model. This problem has been studied earlier 0] by a 
different method which was not based on the SUSY trick. 
The Id Anderson model corresponds to the field theory 
Eqs.QJ,© where the coupling matrix gy = w <5y is diag- 
onal and fly = 6i t j+i + 6i t j—i. To compute the distribu- 
tion function Pi{t) = v-HE) (£ n 5(t - \ip n (i)\ 2 ) S(E n - 
E)) we follow the procedure used in 0,0] for thick quasi- 
ld wires where the sigma-model applies. It starts by the 
exact expression for the q-th moment of this distribution 
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1 



Aw(E) [q- 2)1 Jo 



r+oo 

dss q ~ 2 / dvYi(s,v). 



(10) 

in terms of the function Y(Qi) — JY[e=^i^Q^ e F ^ ■ 
The functional F is given by Eq.QJ, where u> — 
2r\ — > +0 and F kin = Y,i F kin(Qi,Qi+i) = 



-£> 2^J (2- v /StrQ I Q l+1 ). 

As usual 0] in the Id chain of the length N the func- 
tion Yi — WiWn-i e _L ^ i ' is found from the recursive 
relation: 

W j+l (Q) = [ d^{Q')e- F ^^- L{ ^ Wj(Q'), (11) 



where L(Q) = fStrQ 2 + StrAQ - z^StrQ, and the 
boundary condition Wj—\{Q) = 1 is assumed. 

One can show that the integral in Ea. <|llfl does not 
depend on the Grassmann variables x an d on the angle 
ip. Thus Wi and Yi are functions of two (Ai, A2) of the 



seven variables only. This property is based on the form 
of the parameterization Eqs.Q-JHJ and similar to that 
known in the sigma-model approach 0] , where the 
corresponding functions depend on the non-compact Xb 
and the compact Xf angles only. In the limit 77 — > +0, 
Wi remains a function of two new variables v = Ai — A2 
and s = rj(Xi + A2). At this point the difference in the 
symmetries of the target spaces in our approach and the 
sigma-model becomes important. In the latter case Wi 
depends on one rescaled variable ( 2r?) Ab only, while the 
dependence on Xp disappears |3j,|5|,|2lJ. The fact that the 
recursive relation Eq. (|ll(l and the Fokker-Plank equation 
that follows from it depend on two variables is generic to 
all methods in the theory of Id disordered chains 
thus demonstrating that the SUSY approach is the min- 
imal description. We will see that the dependence on 
v reflects statistics of eigenfunction oscillations at scales 
of the order of the de Broglie wavelength hp 1 which are 
essentially different in a thick wire and in a Id chain. 

After integration over the angle if' and over 
the Grassmann variables according to a convention 

/ Xr/A d X' R /A = and g° in g to tnc limit f ~~* onc 
obtains an exact recursive equation: 



W j+1 (s,v) 



s 2 
2^ 



dv' 



ds' 

(Pji 



Wj{a'J) (12) 













X COS 




(7- 




e 













■ v' 2 +iEv' 



Similar equations have been derived earlier |23j, |24j in 
study of the localization transition on the Bethe lattice. 
Ea. (|12|) is valid for any strength of disorder w and for 
an arbitrary energy E. However in the limit of weak dis- 
order w <C 1 it can be drastically simplified using the 
scale separation kp 1 <C £ioc where £; oc is the localiza- 
tion length. To this end we first consider the disorder- 
free case. Setting w = in Ea. (|12|) one notices that 
the Fourier-transform W^s'^q') = J dv' W 3 -,(s', v') e lv ' q ' 
naturally arises. Remarkably, the equation for Wj(s, q) 
turns out to be an algebraic and not an integral equation: 
W j+ i(s, q) = jt e~ sq2 W 3 (sq 2 ,E - |) . This functional 

equation can be simplified if we introduce new variables 
s = z cos 2 (0/2 + k) and q — cos(</>/2)/ cos((/>/2 + k) and 
a new function: 



9j(z,<f>) = Wj {s{z,4>),q{z,4>)) 



sin k 



2cos 2 (</>/2 + fc) 



, (13) 



where k defined by the relation E = 2 cos A: is a momen- 
tum for plane waves in the tight-binding model. 
In new variables we obtain: 



_ e -z cos 2 (0/2) 



$ j (z,4>-2k). (14) 



One can see that the variable (f>/2 is related with the 
phase of the wave function as it changes by k in passing 
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from a lattice site j to the neighboring one. For an ir- 
rational k/w as one moves along the chain at a distance 
iioc ^> F the phase 4> sweeps with the constant density 
through the entire interval [0, 27r]. Therefore <fi is the true 
"fast" variable in contrast to the "slow" variable z. 
Thus one comes to an idea |221 of averaging over (j) and 
replacing <F,-(z,(/>) in Eq. l|13l) by its average (the "fast 

phase" ansatz): <&j(z) = §j(z,<j>) = ^ dcj) $j(z,4>). 

Now we take into account weak disorder by expanding 
e~f v ' 2 = 1 - H v '2 = i + w d 2i dq a in Eq .JJ2l. Wc 

also assume z ~ w <C 1 and expand e _z cos ^/ 2 ) « 
1 — z cos 2 (^/2). Then doing the Fourier transform 
over v and v' , switching to the new variables z,<j) and 
averaging over one obtains for $j+i(z) — &j(z) — 

ciz$j(z) + c 2 z 2 ^p- + c 3 ($j - z^§f}> where a = 

-^P) = -1/2, c 2 = ^(1- cos^) = 
and C3 = -r^a t (1 + cos0)(l — 2cos</>) = 0. Then switch- 
ing to the continuous spacial variable r = j/£i oc with 
&oc = 2s " k and the variable x = z£i oc /2 we obtain the 
following Fokker-Plank equation subject to the boundary 
condition $ T _ (a:) = 1: 

-^=(^-zj* T (*)- (15) 

The fact that C3 = depends crucially on the assump- 
tion on the homogeneous distribution of the fast phase 
4> which is related with the single-parameter scaling. As 
a matter of fact the distribution of phases F(<fi;E/w) 
depends on energy and is homogeneous F(<j>; 00) = ^- 
only for irrational k away from the band center E = 0. 
At the band center F(0;O) oc (3 + cos</>)~ 1/2 is dif- 
ferent [23 and we obtain c\ = 1 — 2E/K w —0.457, 
£2^^ = (-16/3 + 8-E/iT) w 0.494 and c 3 5isiA = 
(-26/3 + VIE jK) w 0.075, where £ = ^(l^ 1 ) and 
ii" = if (1/2) are elliptic integrals. Thus at the band 
center \E\ <C iu Eq.JTSJ should be modified by including 
the combination <I> — xd$/dx proportional to C3. 

Finally we have to switch to new variables in Ea. i|10|) 
which involves the combination Wj(s, q) Wpf-i(s, E — q). 
After replacing W by $ according to Ea. l(T^|l . using the 
fast phase ansatz &j(z, <f) — > $ T (^) and integrating over 
4> one obtains for \E\ ^> w: 

I q (r)= ( 9 2 } " J o dyq(q-l)y«- 2 Y T (Zio C y), (16) 

where we used the relation 2ttv(E) — dk/dE = 
1/(2 sin fc) and the definition Y T (x) — Q T {x) $ n t (x). 

The result given by Eas. l|15|l . (|16|l is in a perfect agree- 
ment with the one obtained in Ref.|l8l| by a completely 
different method. 

Remarkably, the integral in Ea. (|16(l and the Fokker- 
Plank equation Ea. (|15fl coincide completely || with the 



corresponding equations [3J, |5j] determining statistics of 
wave functions for a multi-channel disordered wire. The 
q-dependent factor in front of the integral in Ea. lllfijl is 
due to the difference in statistics of oscillations of almost 
ballistic (in a strictly Id system) and fully chaotic (in a 
thick wire) wave functions at a small scale of the order 
of the de Broglie wavelength. This coincidence is a man- 
ifestation of the single-parameter scaling and it breaks 
down near the band center \E\ -C w. 
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